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Gaussianity is the very base for derivation of the cosmological parameters from the CMB angular 
power spectrum. Non-Gaussian signal, whether originated from experimental error or primordial 
source, could mimic extra power in the power spectrum, thereby leading to a wrong set of parameters. 
In this paper we present a new Gaussianity analysis of the derived CMB singals. It is based on 
the trigonometric moments of phases, which can be concluded with a "mean angle", allowing us 
to see the global trend of non-Gaussianity of the signals. These moments are also closely related 
to the Pearson's random walks. We apply these analyses on the derived CMB maps and their 
different morphologies manifest themselves through these functions. We also comment on rotational 
invariance of the trigonometric moments of phases as a non-Gaussianity test. 

PACS numbers: PACS numbers: 98.80.Hw, 52.34.Mw 



I. INTRODUCTION 

One of the main goals of the ESA Planck mission is to produce the maps of CMB temperature anisotropics and 
polarization and the power spectra, from which the accurate values of the cosmological parameters can be derived such 
as baryonic and dark matter density f2 m , the dark energy density fl^, the optical depth of reionization r r and so on. 
A crucial requirement that enables us to derive the cosmological parameters from the temperature and polarization 
power spectra Ct(£) and C p (£) is that the statistical properties of the primordial CMB signal should be Gaussian. 
Should the primordial CMB signal possess non-Gaussian origin such as quadratic non-linearity in the gravitational 
potential (Komatsu et al. 2003, see also Wandelt et al. 2004, Mattarese et al. 2004), the connection between Ct(£), 
C p (£) and the cosmological parameters needs an additional investigation which seems to be non- trivial. One can 
illustrate the importance of non-Gaussianity of the CMB signal by assuming that at some range of multipoles, say, 
I ~ 200 the ai m coefficients of the spherical harmonics expansion of the anisotropy AT are highly correlated. Without 
comprehensive testing of non-Gaussianity of the map, these correlations can easily mimic the first acoustic peak in 
the Ct(£), leading to a wrong conclusion on the properties of the CMB and the cosmological parameters. Preparation 
and implementation of sensitive non-Gaussianity tests on the anisotropy and polarization maps is therefore pivotal 
for the Planck mission. 

After the release of 1-year results of the Wilkinson Microwave Anisotropy Probe (WMAP ) [3-5, 23, 24] in the papers 
by Chiang et al. (2003), Park (2004), Eriksen et al. (2004), Hansen et al. (2004), Larson and Wandelt (2004), Land 
and Magueijo (2004), Roukema et al. (2004), Schwarz et al. (2004), Chen et al. (2004), various kinds of non-Gaussian 
features are detected in the WMAP derived maps. These structures of the WMAP signal in combination with phase 
analysis of the whole sky maps, derived from the WMAP data sets (Chiang et al. 2003, Naselsky, Doroshkevich 
and Verkhodanov 2004) and wavelet analysis (Vielva et al. 2004, Cruz et al. 2005) clearly show that WMAP signal 
contains non-Gaussian features, which could be primordial origin (Eriksen et al. 2004) or related with foreground 
residues (Naselsky et al. 2004, Chiang and Naselsky 2004, Dineen and Coles 2004). The exercise of finding non- 
Gaussian peculiarities in the WMAP CMB maps is also crucial for the upcoming Planck mission, paving the way 
especially for reconstruction of the signal in the Galactic plane area, in order to obtain the whole CMB sky for both 
temperature anisotropics and polarization. The method of testing non-Gaussianity presented in this paper operate in 
the multipole coefficients ai m and related to phase analysis, developed by Chiang and Coles (2000), Chiang, Naselsky 
and Coles (2004), Coles et al. (2004), Naselsky et al. (2004a,b), Dineen, Graca and Coles (2005). It is applicable for 
both temperature anisotropies and E component of polarization multipole coefficients. 

The purpose of this paper is to show that any correlations of phases for CMB signal determine the morphology 
in the space of phases. We generalize the method of phase correlations by Naselsky, Doroshkevich and Verhodanov 
(2003) for testing of the phase coupling in the CMB maps. Particularly, we introduce the mean angle O^ for each 
multipole £, averaged over all m-modes in order to check whether the distribution of Qg is uniform (as it should be for 
Gaussian signals), or, if not, what possible preferable directions for each miltipole could be. Note that this statistic 
of phases has a natural explanation in terms of Pearson's random walk problem (Pearson 1905) for ^ m ai m and can 
be easily generalized for more complex correlations of phases. The importance of such approach is recently pointed 
out by Stennard and Coles (2004). 

For illustration of sensitivity of the method we use the foreground cleaned map (FCM) and Wiener filtered map 
(WFM) derived by Tegmark, de Oliveira-Costa and Hamilton (2003) (hereafter TOH) foregrounds cleaned map 
(FCM), the ILC map by WMAP science team ( [38]) and the ILC map reproduced by Erikson et al. (2004) (hereafter 
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EILC map). All these maps contain some features of the foreground residues, non-uniformity of the noise, Galactic 
plane substraction . . . etc. To avoid confusion, below we call these maps the CMB maps, which represent different 
morphology of the primordial (Gaussian) CMB signals with some of the non-Gaussian features. The main task is 
to show how any detected non-Gaussian features relate to known properties of the non-Gaussian components of the 
signal. All these maps include Galactic plane contamination, which is excluded by different kind of masks in order 
for estimation of the power spectrum of CMB anisotropies. 



II. PHASE STATISTICS AND THEIR APPLICATION FOR THE CMB ANISOTROPY PROBLEM 

In this section we formulate the problem: how do the anisotropies of the CMB signal manifest themselves in the 
phase correlations? To answer this question, we recap some of the basics of Gaussian random field. The CMB 
temperature fluctuations on a sphere can be expressed as a sum over spherical harmonics: 

oo I 

AT(6,v)=J2 E Wi,n\e^Y hn (e^) 7 (1) 
e=o m =-e 

where \ai m \ and 4>i m are the moduli and phases of the coefficients of the expansion. 

In practice, we use the Healpix package [21] to decompose each of the following derived CMB maps (ILC, FCM, 
WFM, EILC) for the coefficients of spherical harmonics at m , extracting the phases for analysis. 

A homogeneous and isotropic CMB Gaussian random field (GRF), as a result of the simplest inflation paradigm, 
possess Legendre polynomial expansion modes whose real and imaginary parts are Gaussian and mutually independent 
[6, 7]. The statistical properties of such a field are then completely specified by its angular power spectrum C% mb , 

{aj™ b (a cmb )} m ,) = Cr b S ee '6 mm> . (2) 
In other words, the Central Limit Theorem guarantees that the field is Gaussian if their phases 
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are randomly and uniformly distributed at the range [0, 2ir]. The method we propose is therefore based on the so- 
called "random phase hypothesis" , our null hypothesis for a Gaussian random field, by testing any phase correlations 
between modes. For this purpose we apply circular statistics [20] on phases of the spherical harmonic coefficients 
(Naselsky et al. 2003; 2004). We examine the cross-correlation of phases between modes {£, m) and (£ + At, m + Am) 
for all £, m values. The cross-correlation of phases vanishes by definition for pure Gaussian CMB signal, whlist for 
non-Gaussian components it should display non-trivial significance. Note that we do not examine any specific kind 
of non-Gaussian signals. In reality, all different sorts of non-Gaussian components can propagate into the data, such 
as systematic errors and foreground residues, in addition to primordial non-Gaussianity (if exists). 

The basic idea is to introduce trigonometric moments of phases, which minimize the contribution of the non- 
correlated Gaussian tail and maximize the non-Gaussian tail of the phases (Naselsky, Doroshkevich and Verkhodanov 
2003; Naselsky et al. 2004). To simplify the analysis let us describe the phase correlations in £ — m plane separately, 
choosing orthogonal directions At ^ 0, A m = and At — 0, A m ^ and throughout this paper we demonstrate 
phase coupling in these two directions. Following [29, 30] we define the following trigonometric moments for the 
direction At ^ 0, A m = 0: 

£ 

C(£, At) = cos (^£+Ai, m - ®e, m ) ; 

m— 1 

t 

S(£, A e ) = sin {^t+At,m - $e,m) ; 

(4) 

and for the direction At = 0, A m ^ 

^max 

C(m, A m ) = Y cos _ ®e,m+Am) ; 

l=m 
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S(m, A TO ) = sin (* <>m - $£, m _Am) ; 

(5) 

where f max is the maximum ^ putting into calculation, and can be the phases of the spherical harmonics 
of two maps for cross correlation, or those of the same map (\IV TO = &e m ) for auto-correlation. The concept behind 
such statistics is clear. If the phases of the signal are highly correlated (^e,m — ► &e+Ae,m), for the Sng statistics we 
get S — > 0, while C — > 1. We therefore can further device the "mean angle" along m and £ direction as 

e <-'' A ™> ='»"-' Ira < 6 > 

respectively, both defined in [0, 2tt]. Thus, for the pure Gaussian signals the mean angles should be uniformly random 
along either the £ (after summation over m) or m (over £) direction. Non-Gaussian signals in the CMB data manifest 
themselves as asymmetry with respect to ir or non-random. Moreover, for non-Gaussian signal we expect to find 
cross-correlation between Q(£, Ag) and Q(£ , A^), and between 6(m, A m ) and 6(m , A m ) for different £, £ and m, m 
modes. 



III. TRIGONOMETRIC STATISTICS AS AN EXAMPLE OF PEARSON'S RANDOM WALK 

One of the interesting properties of the trigonometric statistics is related to the idea suggested in Stannard and 
Coles (2004), in which it is pointed out that the statistics of the ag m coefficients in the form of = J2 m >o aim ^ s 
essentially the same as those of "Rayleigh flight" , whose properties are well developed. However, for the trigonometric 
statistics proposed in the previous Section, we need to apply Pearson's random walk statistics, taking into account 
the investigation of the random walk with given (with fixed length of steps) by Pearson in his famous letter to the 
Journal Nature [32]. Following Hughes (1995) and Stannard and Coles (2004), we introduce the random variable 

m=M m=M 

Gi (M)= £ y^= £ (7) 

m=l 1 1 m=l 

the real and imaginary parts of which are essentially following the properties of the Pearson's Walk in the complex 
plane after M steps. Note: that the walk can be in £ direction for G m (M). For all the vectors G(M) the starting point 
is zero and the length of each step is unity: \ae m \ — 1- As one can see from Eq.(7), vector Gg(M) has a projections 
on x and y axes ^cos\I^ m and ^sin\I>£ m , which are the trigonometric moments introduced above. The probability 
density function (PDF) for the position of the walks after M — £ steps satisfies the equation (see Hughes (1995) for 
references) 



Pm+i(v) = J p M +i(r - r )P M (r )d 2 r (8) 

where PM( r ) is the probability density function that the vector G{M) lies in an infinitesimal area centered on r on 
the M-th step. For isotropic two-dimensional random walk, in which each step has the same length a=lwe get for 
p M (r) (Kluyver 1906) 

PM=i(r) = ^S(\r\-l) (9) 

and the corresponding PDF for M steps is 

1 f°° 

p M(r) = —J^ dv vJMA) IM")} M , (10) 

if M > 2 and Pi(r) = pi(r) for M = 1. If M > 1 the PDF has Gaussian asymptotic (see the references in Hughes 
(1995)) 
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So, for M ^> 1 the corresponding variance of the random process is a 2 = M/^/2. We would like to point out that this 
result can be obtained directly from the circular variables C and S by the following way. Let introduce the random 
variable R 2 M — C 2 M + S 2 ^ (Fisher 1993). From Eq.(5) the upper limit now is M and = 0, one can have 



M M 



Rm=Y, C0S (*^ - ) = M + S C0S (*^™ - )' ( 12 ) 

For ^> 1 and non-correlated phases the last term in Eq.(12) is in order of M _1 and R 2 M = M. 

If the vector Gm{£) defines the distance from zero point, at what direction do the walkers move in the plane after 
M steps? The answer to this question can be obtained from circular statistics, namely, by Q(£, A^). Moreover, one 
can formulate the problem of two walks, which randomly move on the plane from the same initial point C and after 
M steps come to some points A and B. What is the probability to find both walkers at the vicinity of that points, if 
the distance between A and B is r? More importantly, if the steps of two walkers are correlated, what is the distance 
between them after M steps? The list of problems can be significantly extended, but using Pearson's walk model we 
can reformulate some quantities of the trigonometric statistics in terms of random walk approach. However, some of 
the statistical estimators can be easily derived using trigonometric moments, which we will demonstrate in the next 
section. At the end of this Section we would like to point out that the properties of the 8(m, A m ) statistics can be 
described in term of random walk model, but now the "walker" made the steps in I direction starting from m = € max 
and down to m — 1 having corresponding variance R 2 n — £ max — m + 1. That sort of statistics allows us to investigate 
the morphology of the space of phases for each fixed value of m in direction m < I < £ max . 



IV. TRIGONOMETRIC STATISTICS IN APPLICATION TO WMAP DATA 

In this section we will investigate the asymmetry of different CMB maps derived from the WMAP data sets by 
Bennet et al. (2003) (ILC map), the Foreground Cleaned Map (FCM) and Wiener Filtered Map (WFM) by Tegmark 
ct al. (2003), the Internal Linear Combination map reproduced by Erikson et al. (2004) (EILC). All these maps are 
derived by different methods of CMB and foreground signal separation, applying different parts of the sky with masks 
and disjoint regions and they have different angular scales for the CMB image reconstruction. In the analysis which 
follows, we are not going to discuss the limitation of the methods applied by the authors of ILC, FCM, WFM and 
EILC, or use these maps as possible illustration of different statistical properties of the signal. Our goal is to show 
the sensitivity the trigonometric and Pearson's random walk statistics in application to different morphologies of the 
maps. For simplicity we will continue to use the term CMB map, but the reader should bear in mind that we use 
these maps for their illustrative character. 

For example, the FCM derived with some smoothing by the Gaussian filter of the signal at angular scale 8 < 1°, 
where the amplitudes (but not the phases!) of the signal is vanishing within that scale. The WFM has some additional 
linear (Wiener) filter which was applied separately for each disjoint region of the FCM using the WMAP best-fit power 
spectrum Ci. The question is: how can we compare these maps, which have different resolutions and different error 
bar level after foreground and noise separation? 

One may argue that the best way to compare two maps is simply taking the difference between 2 maps (pixel by 
pixel). However, due to different resolutions and different errors, taking the difference of 2 maps is by no means a 
good estimator. To elucidate this point, we show in Fig.l the function w{tj = W?(t), where W{£) is the window 
function for all WMAP K-W bands. The properties of the w{i) . In each band there are two parts of singal, Si and 
S2, where Si is the sky signal So convolved with the beam B: Si = So ® B, and S2 is the noise. The phases of the 
ai m harmonic coefficients are therefore 

, _ s l<5i^|sin* £ 5l m + |S 2 , !m |sin*f m 

tan*?.,, = n g— (13) 

' |S 1/>m | COS *f m + 152/^10)8*^ 



If |Si^ m | ^> \S2t, m \> which is a good approximation for the low I range, the properties of the phase *f TO are the 
dominated by signal Si. However, when |Si^ m | -C m | (for the high I range), the phases of the signal *f m are 
determined by the phases of the noise. Using these asymptotics, we can conclude that the range of multipoles for 
which \Si£ m \ ~ |S 2 ^ m | is the range of transition (£ ~ £*) from the phases of the signal from the sky to the phases 
of the noise (see Chiang ct al. 2002 for details). With different beams and window functions, the transition range is 
different for different bands. 

We propose another way for comparison of the different CMB and non-CMB maps, based on asymmetry parameters. 
Namely, for each map, using coefficients of the spherical harmonics expansion, we derive the phases §\ for comparison 
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FIG. 1: The window function w(£) and 1 — w(£) functions for the K-W bands versus I. The solid lines and dash-dotted lines 
are w(£) and 1 — w(£), respectively, from top to bottom for WMAP K, Ka, Q, V and W bands. 

using auto- and cross-correlation of the phases. Below we present some of the results of the analysis. We include in 
our analysis not only phases of the derived CMB maps, but phases of the K-W bands of the WMAP data sets in 
order to compare the properties of the CMB and original data. In Fig. 2 we plot the trigonometric moments for the 
phases of the WMAP Ka-W band signals. These are supposedly pronounced non-Gaussian signals due to the Galactic 
contaminations, and Q a (£, A^), 9 a (™, A m ) display the highly concentrated regions, in particular for 5g = 2, therefore 
indeed non-Gaussian. 



V. PHASE AUTO-CORRELATIONS FOR THE DERIVED CMB MAPS 

In this section we discuss the properties of the CMB maps derived from the WMAP data sets by combination of 
the K-W bands maps. As is mentioned in Introduction we use the ILC, FCM, WFM and EILC and the corresponding 
ag m coefficients for all range of multipoles (^ max = 512). However, none of these maps represent pure CMB signal up 
to this range of I because of different restrictions of the methods applied for separation of the CMB signal, noise and 
foregrounds. Meanwhile, these maps allows to test the properties of the trigonometric statistics of the phases, which 
as we believe can be useful for testing the upcoming WMAP data and the Planck data. 

The statistics we use are Q a (i, Ag), 9(m, A m ) for A^ = 1, 2 and A m = 1,2. 

For the random signal all these functions arc uniformly distributed by definition at the range [0, 2tt], while for the 
non-Gaussian signal they should display clusterization of the phases. To compare between the phases of maps, in 
Fig.3 we plot e(^A^), 9(m, A m ) functions from left to right ILC, EILC, FCM, WFM, respectively. We also would 
like to draw attention on the morphology of the Q a (£,Ag), 6(m, A m ) functions for Ag = 1,2 and A m = 1,2. The 
correlations for the WFM are discussed in Naselsky et al. 2004 and are independently found by Prunet et al. (2004) for 
the WMAP signal. However, for the ILC map the 6(m, A m ) statistics have displayed some non-Gaussian character, 
where for A m = 1, 6 are more clustered between n and 27T and for A m = 2, it is more concentrated around n. 

VI. CROSS-CORRELATIONS BETWEEN THE DERIVED CMB MAPS 

In this Section we would like to demonstrate how <d c (£, Ag), 6 c (m, A m ) statistics can be used for comparison of 
different CMB maps derived by different methods, including the linear (Gaussian or Wiener) filters. These linear filters 
significantly destroy some information at specific ranges of multipoles, but not the phases of the signals. Convolving 
with a Gaussian filter with characteristic scale £ = £ G , all the amplitudes ag m for £ » £q become vanishing while all 
the phases preserve all information about the morphology of the CMB map. That is one of the reasons why phase 
analysis is essential for testing of the statistical properties of the derived CMB maps, as it clearly illustrates different 
outcomes of the cleaning methods used. To compare the phases of ILC, EILC, FCM and WFM we implement Eq.(5), 
in which the phase 'F now corresponds to one map (for example, ILC) and the phase $ corresponds to another (for 
example, FCM). We exploit the simplest fact that if the phases are the same, in terms of Q c (£,Ag), 9 c (m, A m ) 
statistics, they behave like strong non-Gaussian signal. In Fig. 4 we plot the cross-correlators <d c (£,Ag), 9 c (m, A m ) 
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FIG. 2: Auto-correlation of phases in terms of Q a {i, At), a (m, A m ) for the WMAP (from left to right columns) Ka, Q, V 
and W bands. From top to bottom are A e — 1 and A m = 1, Ai — 2 and A m = 2, respectively. 



by pairing CMB maps: from left to right are between ILC and EILC, ILC and FCM, ILC and WFM, and between 
EILC and FCM. 

Comparing the phases of ILC and EILC maps one can see significant ±7r anticorrelations for I > 400 at A£ = 0, 
while for other pairs they are negligible. And at Am = all pairs start deviating from total correlation for m > 300, 
indicating that the main difference in morphology of these derived maps starts from the multipole range £ > 400. 
More interestingly, for the ILC and EILC maps the phase correlations have systematic shift by ±7r. Note that while in 
terms of C (£, Ai = 0) and 6 c (m, A m = 0) statistics the signals are nearly the same for I < 400, they have different 
auto-correlation (see Fig. 3), which is why the Q c {£, Ag — 1, 2) and 6 c (m, A m = 1, 2) statistics are different. 



VII. DEVIATION FROM RANDOMNESS IN TERMS OF PEARSON'S WALK STATISTICS 

As is mentioned in Section III, the properties of the trigonometric moments can be understood in terms of the 
Pearson's random walk. It provides some additional information about auto and cross-correlation of phases even if 
the deviation after M steps from the estimated Rm is 2 — 3 times bigger than for the Gaussian estimator Rm = M 1 / 2 . 
Although the 6 C (£, A^) and c (m, A m ) functions show the gross behavior for a certain separation, the random walk 
of phases can display how the phases correlate step by step. 

In Fig. 5, we present the random walk for Gaussian random signal, ILC, FCM and EILC with fixed separation 
AC. = 2 for £ = 350 and Am = 2 for m = 6. We choose the particular case of A£ = 2 for I = 350 because the resultant 
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FIG. 3: Auto-correlation of phases in terms of © a (£, A^), a (»n, A m ) functions from left to right columns for the ILC, EILC, 
FCM and WFM, respectively. From top to bottom rows are Ae = 1, A m = 1, Ag, = 2 and A m = 2, respectively. 

angles for EILC and FCM are close to n, as shown in Fig. 3 the 3rd row, the same reason for Am = 2 for m = 6 
for ILC and FCM. The Gaussian a values are calculated when we set the L max = 350 and M max = 350. The walks 
projected in two axes can be written as 

M 

C<?=350,A*=2(M) = c Os(*£ + A £ ,m ~ ^ Lm) 

m—1 
M 

Sfc350,A f = 2 (M) = ^ sin(^ +A „m ~ */, m ) (14) 
m—1 

and 

L 

C m=6 ,A ro =2(i) = COS (*<>"» - **,m-A m ) 

L 

S m=6iAm = 2 (i) = J! Sin (^,m - */,m-A ra )- (15) 

As one can see from Fig. 5 the walks of the derived CMB maps are intrinsically different from the random signal. 
For ILC, the walk with the phase difference Am = 2 reaches 3cr, so do the walks with A£ = 2 for EILC, FCM and 



8 



100 mo am wo am 




im a an n:- ra 
tmffilpclJB [ 



af 

* E 

* 1 
£ at 



III ffi* .3)} IKi 



I 9 ■ 
4) 

ITS 



« 



1!« !D) J(0 #! IB] 




1(M SKl 300 U0 Ml 



=1 r 
■ 5 
* 4 

I. 

V B 

1 






ido zn am 400 ehu 




1D0 MO 5] [1 400- HO 
Bqpatilbil number -» 




-1 



El 



i 

-'ii 



1 aniiffn 



1H> Zfl JZE 1 -KD S¥J 



in sjd -HD nw 



, E 




idd BOD SB- 40D HB 
iquulUliiaibir ■» 





1W SM W -WJ PW 



IM 2W WH *W WC- 



1WHW 3TO -Wf- WJ 
nrUHjmlri mMLtv ' 



1OT .3TO -HM -W BW 




in ani' m *oa mi 




1. 

2„ 




id ma «] an 




FIG. 4: Cross correlation of phases in terms of Q c {(.,Ae), c (m, A m ) functions between a pair of maps. From left to right 
columns are between ILC and EILC, ILC and FCM, ILC and WFM, and between FCM and EILC. From top to bottom rows 
are At = 0, A m = 0, A* = 1, A m = 1, At = 2 and A m = 2, respectively. 
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FIG. 5: Random walk of ^ M exp(iA3>). In each panel the x and y axes are C(M) and S(M) (top row), and C(L) and 
S(L) (bottom), respectively and the circles (starting from inner one) represent l,2,3,4er . . . of Gaussian confidence levels for 
maximum M max = 350 and Z/ max = 350. The top row is for Ae = 2 of ^ = 350 and bottom for A m = 2 of m = 6 with 
imax = 350. From left to right columns are Gaussian random signal, ILC, EILC, FCM and WFM, respectively. 



WFM. Although the statistics from trigonometric moments cannot provide any information about non-Gaussianity 
for a single t or to, through Pearson's random walk we can see significant correlation of the phases in any chosen 
direction. The corresponding values of the Rm parameter for the maps are around or above 3<r. 



VIII. ROTATIONAL INVARIANCE OF THE TRIGONOMETRIC MOMENTS OF PHASES 

The previous analysis of the trigonometric moments as estimators of non-Gaussianity of the anisotropy and po- 
larization maps dealing with phases of ag m coefficients consequently depends on the reference system of coordinates. 
Obviously, for different coordinate systems these ae m coefficients and the corresponding phases are different. The issue 
which we would like to discuss in this section is how phase correlations depend on the reference system of coordinates 
(M. Hobson, private communication, see also Coles et al. 2003). We concentrate on the following question: are the 
phases correlators, in particular, the trigonometric moments we introduce, rotationally invariant? And if not, how 
significantly such non-invariance can transform the conclusion about non-Gaussianity of the maps? To answer these 
questions we need to know how to transform a given set of at m of a given coordinate system to new coefficients bi mi 
which corresponds to the new coordinate system rotated by the Euler angles a, /3, 7. Following the general method 
(Varshalovich, Moskalev and Khcrsonskii 1988) and taking into account properties of the spherical harmonics (see 
Coles et al. 2003) we get 

m 

where D l , (a, (3, 7) is the spherical harmonic decomposition of the Wigner function D(a, /3, 7). The coefficients 
D e m , (a, (3, 7) should preserve the moduli J2 m \be,m\ 2 = J2 m \ a e,m\ 2 under transformation. This leads to the following 
equations for D l m m , (a, (3, 7) 

E i^»i 2 -EE C'( a '^^^"( a ^^)v a V = E k™i 2 ( 17 ) 

m m m ,m" m 

and 

E D L >m '( a ^^) Da m,m" («,/?, 7) =5 m . im . (18) 

m 

where 8 ' •• is the Kroneker 8 symbol. Without loss of generality we assume that for the initial reference system 
of coordinate ai ;m — exp(i$£ !m ) where $>i ;m is the phase for a given £ and \m\ < L To obtain the trigonometric 
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moments for the system after rotation we define a matrix of correlations: 

G t/ =E 6 ^',™ = E E D L,m'^^) D im" (M,7Kw fl /,.- ( 19 ) 

where £ = £ + A^ and &£ m = |^, m | exp(i^^ iTO ), and \f^, m ) is the phase of the £,m harmonics after rotation. For 
the statistical ensemble of realizations of the same Gaussian random process, average over realizations leads to 
orthogonality of the a^ m coefficients: (a lm >a*, m „) — 5 m > m -S le '. Thus, the correlation matrix has the form G ee > — 

S e gi which is typical for a Gaussian process. Unfortunately, in CMB studies we are dealing with one single realization 
of the sky, therefore such asymptotics can not be achieved. In a statistical sense, however, non-correlated phases 
in a Gaussian signal preserve randomness of the trigonometric moments and could play a role for estimation of any 
non-Gaussianity of the map. 

Let's discuss the properties of trigonometric moments for rotated maps. From Eq.(7) we obtain 

e(£,a,/3, 7 )=tan- 1 — ! '- ■ L (20) 

^{E m E m ', m » ^ m ,(a,/3, 7 )£>^ m .(a,/3,7)exp[i(^ v -^, m -)]} 

Taking into account (see Varshalovich, Moskalenko and Khersonskii, 1988) 

D L, m ' = e M-ima)d e m m , (0) exp(-«m' 7 ); (21) 

where 

i min(<-m,<-m ) f / . _ . 

d l ,{0)= {£ + m)\{£-m)\{£ + m)\{£-m)\ 2 V {-if ~ m - k - ,"V 

y ' {£-m -k)\{k + m + m)\ 

X (COS 0f k + m +™ (gin pft-U-™-™' ( 2 2) 



A- 



H 


E ro ', m -Wi( m .C9)sin($, v 


- $ £',m" + 7 (™' 


-m")) j 




E m ', m '-^; m ..(/3)cos($, v 







we obtain the following formula for 6(€, a, /?, 7 ): 

cvir , „ 

e^,a,/3, 7 )=tan- 1 --4-— — - : — ^ (23) 

where 

^0^=E<w(/3)<^(/3) (24) 

m 

As one can see from Eq.(23) the @(£, a, 0, 7 ) statistic is rotationally invariant if /3 = and 7 = 0. It is practically 
rotationally invariant even when = but 7 7^ 0. For that case as Eq.(23) the main contribution to the summation 
in both the nominator and denominator corresponds to m — m and Q(£,a,0 = 0, 7 ) ~ Q(£,a = 0, = 0, 7 = 0), 
where ©(^, a = 0, /? = 0, 7 = 0) corresponds to the initial reference system. However, if 7^ 0, the rotational 
invariance is broken and the statistical properties of Q(£, a, 0, 7 ) and &(£, a — 0, — 0, 7 = 0) are different for given 
multipoles m. To investigate how significantly the &(£, a, 0, 7) statistics can be changed we resort to the Monte 
Carlo simulation presented below. We use the GLESP code (Doroshkevich et al. 2003), which provides map rotation 
for arbitrary angles 9 and (f>. In Fig. 6 we show the rotation of the FCM map by 9 — 55° and <j> = 40°, and calculate 
the trigonometric moments shown in Fig. 7. 

As one can see from Fig. 7, the non-Gaussianity of the FCM map is preserved in different form in terms of the 
mean angle statistics @(£, a, 0, 7). The residuals from the Galactic plane contributes significantly to Q(£, a = 0,0 = 
0, 7 = 0) statistics, starting from high multipoles (£ > 200). After rotation, the Galactic plane is shifted causing more 
correlations for Am = 1. In fact, non- invariance of phases provides a new way to test any residues, including from the 
Galactic plane, by rotation of the map by 9 = n/2, for which any alignment in the structure in the original orientation 
is now perpendicular to the previous one. Once again we would like to emphasize that by definition for pure Gaussian 
signal any rotations should not lead to correlation of phases, which can be used for estimation of significance level of 
the non-Gaussianity. 
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FIG. 6: Rotation of the maps. The left is the FCM map and the right is its rotated map by 6 = 55° and <f> = 40°. 
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FIG. 7: The mean angle statistics for auto correlations for the rotated FCM map as shown in Fig. 6. From left to right are 
At — 1, A m = 1, At = 2 and A m = 2, respectively. These are to be compared with the 3rd column in Fig. 3. 



IX. CONCLUSIONS 



In this paper we have presented a new method for analysis of non-Gaussianity for CMB maps, is further development 
from the trigonometric moments of phases of the spherical harmonic coefficients ag m can be concluded to form mean 
angles that enables us to see the global behavior of the "randomness" of the phases. Furthermore, the trigonometric 
moments are closely related to the Pearson's random walks and the mean angle is the resultant angle of the 2 
dimensional random walk on the complex plane. The properties of such random walks are well-developed and they 
are used to display the phase correlations step by step. We apply these analyses on the derived CMB maps from 1-year 
WMAP data. These CMB maps are derived from different foreground cleaning methods, hence their morphologies 
are somewhat different. As phases are closely related to morphology, our analyses on phases not only demonstrate 
the existence of non-Gaussian residuals among these CMB maps, but also reveal the different morphologies of these 
maps. With the upcoming ESA Planck mission, CMB maps with higher resolution and sensitivity will be derived by 
different foreground cleaning methods, for which testing the Gaussianity will be imminent. The methods and analyses 
we present in this paper for the WMAP data can be applied straightaway to new datasets by the Planck mission. 
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